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Breaking anchored droplets in a microfluidic Hele-Shaw cell 
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We study microfluidic self digitization in Hele-Shaw cells using pancake droplets anchored to 
surface tension traps. We show that above a critical flow rate, large anchored droplets break up to 
form two daughter droplets, one of which remains in the anchor. Below the critical flow velocity 
for breakup the shape of the anchored drop is given by an elastica equation that depends on the 
capillary number of the outer fluid. As the velocity crosses the critical value, the equation stops 
admitting a solution that satisfies the boundary conditions; the drop breaks up in spite of the neck 
still having finite width. A similar breaking event also takes place between the holes of an array of 
anchors, which we use to produce a 2D array of stationary drops in situ. 
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One of the underlying drivers in microfluidics is to 
miniaturise the standard multiwell plate and transform it 
into an integrated programmable device, while replacing 
the hundreds of wells typical today with thousands or 
more nanoliter-scale compartments. Early success was 
achieved by using deformable chambers that could be 
opened or closed using an external pressure source [T] . In 
parallel droplet-based microfluidics continues to attract 
ever increasing interest since it provides an elegant way 
to encapsulate an initial sample into a large number of 
independent micro-compartments. However, the tradi¬ 
tional droplet production and manipulation methods all 
rely on the drops flowing in a row in a linear microfluidic 
channel. Studying the contents of these drops mm is 
therefore more akin to flow cytometry than to multiwell 
plates. 

Recently several groups have shown how to array 
droplets in micro-fabricated traps, particularly in a wide 
two-dimensional ( 2 D) region [HIS]. These devices typ¬ 
ically allow a droplet density of several hundred per 
cm^, far beyond what is currently possible outside mi¬ 
crofluidics. These devices must still be coupled neverthe¬ 
less to a traditional drop production device before these 
are brought to the observation chamber. While the un¬ 
derlying physical mechanisms for these drop production 
devices is now well understood 0171 , they are poorly 
adapted in practice to making a limited number of sta¬ 
tionary drops. For this reason, quasi-two dimensional 
devices have been designed to break an initially large 
drop into stationary sub-droplets that are held in pock¬ 
ets on the side of a sinuous channel [BHIQ], while truly 
2D devices would allow much higher density of trapped 
droplets m- 

In this letter we describe the ability to break droplets 
in-situ in a wide chamber, by pushing them over a truly 
two dimensional array of micro-fabricated traps [ 12 ]. 
We elucidate the physical mechanisms first on a single 
droplet and show that it is well described by a set of uni¬ 
versal curves given by the elastica equation. The drop 


then breaks through a singularity in the curves beyond a 
critical deformation, leading to a well characterized and 
robust size. In addition to its applications for droplet 
arrays, this new route to breaking a liquid interface pro¬ 
vides fundamental insight into the evolution of drops and 
bubbles in all confined geometries, where the traditional 
Rayleigh-Plateau instability is not active. 

The experimental setup consists of wide chamber 
(width w = 2.5 mm) in which the drops can be an¬ 
chored in a central trap, as shown in Fig. l^-b. The 
chamber height h and the anchor diameter d were in the 
range 35 — 50 /im and 75 — 200 /rm, respectively. The 
height of the anchor e is chosen to be of the order of the 
height of the channel. All devices were made out of poly- 
dimethylsiloxane (PDMS, Dow-Gorning). The droplets 
studied here were produced by a flow-focusing junction 
upstream of the chamber (Fig. Ef), designed such that 
the drops were forced to adopt a pancake shape, of ra¬ 
dius R h, in the chamber. Droplets of fluorinated oil 
were produced in a glycerol/water mixture with 2 % SDS 
as a surfactant. The oil viscosities /ioii ranged between 
1.2 and 24 cP, whereas the outer phase viscosity /i was 
varied between 0.89 and 3.3 cP by varying the water to 
glycerol ratio. The interfacial tension between the drop 
and the outer liquid had a typical value 7 ^ 17 mN/m. 

During a typical experiment, a droplet is anchored on 
a surface energy trap and the outer velocity is increased 
step-by-step. A typical measurement of the drop evolu¬ 
tion is shown on Fig. [^, where the smallest width of the 
drop neck is reported. We observe that the neck reaches 
a stationary value when the outer fluid is flowing below 
the critical velocity. At the critical velocity, the neck be¬ 
gins to decrease from a finite value, slowly at first but 
then this decrease is accelerated until the neck width be¬ 
comes equal to the channel height. Because of this loss of 
confinement locally, the Rayleigh-Plat can instability be¬ 
comes active and the neck breaks, as observed previously 
for different droplet production situations 0I71II3]. 

For a given setup and droplet volume, we find that 
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FIG. 1. (a) Schematic of an anchor of diameter d and 

height e, trapping a drop of radius R (no external flow), the 
chamber height is h. (b) Sketch of the microfluidic device: 
droplets are produced with a flow focuser and trapped on an 
anchor centered in the observation chamber, ic = 2.5 mm. 
(c) Typical evolution of the droplet neck width (see inset) 
as the outer flow velocity U is increased. When the crit¬ 
ical velocity is passed the neck never reaches a stationary 
state; instead it decreases until it matches the channel height 
{h = 49 //m), where it becomes unstable by Rayleigh-Plateau 
instability, (d) For a given microchannel, a trapped droplet ei¬ 
ther escapes from its anchor at a critical value of the capillary 
number (greyed area, R < 600 //m) or breaks on the anchor 
(R > 600 /xm), leaving a smaller droplet behind. Experimen¬ 
tal points for droplets of FC-40 in water/glycerol mixtures of 
different viscosities: = 0.93 cP (dark grey), /x = 1.6 cP 

(medium grey), and fi — 3.3 cP (black). SDS was used as a 
surfactant at a concentration of 2% in all cases. 


there exists a critical capillary number Ca"*" = nU'^ 
above which the droplet either escapes from the an¬ 
chor [12] (Fig.j^dl), or breaks on the anchor (Fig.[^d2). 


In this second regime, the drop leaves behind a daughter 
droplet small enough to remain anchored at the breakup 
velocity the daughter droplet lies in the trapped re¬ 
gion of the phase diagram, as indicated by the diamond 
in Fig. Ef- In the rest of this letter we focus on the 
breakup regime exclusively. 
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FIG. 2. (a) Evolution of the htting parameter if with the 

value of Coilh? ^ for different channel heights, trap diameters 
and droplet fluid viscosities. The cloud of data points con- 
hrms that ~ . (b) Photograph of an anchored 

droplet under flow and system of coordinates in which the 
interface r(s) is derived, s is the arclength along the inter¬ 
face and r'(s) is its tangent, (c) Erom left to right : droplet 
shape for increasing outer fluid velocities. Pink dotted lines 
are the best htting elastica shape, (d) Three droplets of dif¬ 
ferent viscosities (/xi = 1.2,/X 2 = 4.1,/xs = 24 cP) but similar 
volumes, for the same value of /x^Z/y/i^: the droplet shape is 
independent of the droplet viscosity. 


The droplet stationary shapes are imposed by the pres¬ 
sure difference between the two sides of the interface. The 
flat microfluidic chamber can be modeled as a Hele-Shaw 
cell so we adopt a two-dimensional depth averaged for¬ 
malism to describe the system. Then, the pressure drop 
in the flow direction x is related to the average flow ve¬ 
locity U according to 

Po{x) = -24^x + Cst, (1) 

where Po{x) is the pressure in the outer, aqueous phase. 
The pressure pi inside the oil drop is constant [nun] and 
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is related to the outside pressure Po{x) by the Laplace 
relation: 


TT 

Pi-Po{x) ='f{K±+ (2) 

where and are the curvatures in the perpendicu¬ 
lar and the parallel planes respectively m- Away from 
the anchor, the curvature in the perpendicular plane is 
assumed constant = 2/h. 

The droplet 2D shape r(5), where s denotes the in¬ 
terface arc-length, is described in Fig. It- The tan¬ 
gent to the droplet interface is obtained by differenti¬ 
ating the position with respect to the arc-length r'(s) = 
(cos^( 5 ), sin^(s)), and the in-plane curvature is given by 
/^il = 0'{s). Differentiating Eq. Q with respect to s and 
using Eq. 0 in the particular case where x = r(s) - ex 
yields an elastica equation: 

d2 

^ COS OCs) = 0, (3a) 

where £ = hy^Ca~^/^ denotes the “visco-capillary 
length” of the problem and the radius R of the unde¬ 
formed droplet is used to non-dimensionalize the problem 
(non-dimensional variables, such as 5 , are denoted with 
a bar). This is a pendant drop equation [T2j[T6]. 

Note that the droplet shape does not depend on its vis¬ 
cosity: the only viscosity entering the problem (through 
£) is the external fluid viscosity. To qualitatively check 
this prediction, we show on Eig. 01 the shapes of three 
different droplets of different viscosities (between 0.77 to 
24 cP) but similar volumes, and at the same value of 
Ca. The drops indeed have similar shapes that cannot 
be distinguished from each other. 

Secondary variables are now defined to ease the inte¬ 
gration of the problem. The cumulative dimensionless 
area swept by the droplet interface, a(s), is defined as 
the solution of 


a^(s) = y(s)x'(s). (3b) 

The drop shape is obtained by integrating Eqs. 0a & 
b with a shooting method. Eor given values of d and Z, 
we use the curvature kq = = 0) as the sole shooting 

parameter and search for the values of kq that satisfy 
the geometric condition y = d ai a = n/2. The droplet 
shape r(s) is then fully defined by the triplet {d, Z, 77o}. 

The calculations of 77o are first performed for different 
values of Z, while keeping d constant, and the whole pro¬ 
cess is repeated for different values of d. This leads to 
a catalog of shapes that can be used to fit the experi¬ 
mental droplet shapes. Eor each experimental condition, 
the best fitting shape is found and its associated triplet 
is called {d/, 77o/}. Since all numerical shapes are ob¬ 

tained using dimensionless variables, it is also necessary 
to re-dimensionalize them and we call Rf the radius of 
the best fitting shape for each fit. There is a very good 


agreement between numerical and experimental shapes, 
as shown on Eig. 0- We also find excellent agreement 
between Rf and i?, with a difference of at most 5% be¬ 
tween both values, the best fits being obtained for low 
values of £. 

A quantitative comparison between theory and exper¬ 
iments is obtained by comparmg the evolution of the 
fitting parameter £j‘^ = R‘fl£‘f with the experimental 
value of Ga^jlr?. The resulting values are plotted in 
Fig.0, which displays the data for different channel ge¬ 
ometries, droplet and outer fluid viscosities, and trap di¬ 
ameters. The data verify the scaling predicted by theory 
^ Ca/d^, with a prefactor (^ 60) that differs from 
the prediction (OG/tt — 31), most likely because of dy¬ 
namic surfactant effects [12]. Taken together, the above 
results show that the droplet shape is well described by 
the elastica for forcing values below the breaking thresh¬ 
old. 
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FIG. 3. (a) Branches of solutions of the elastica shown in 
the plane each corresponding to a given value of 

d. The stable part (reps, unstable) is denoted with a black 
line (resp. gray line). For a given d, the critical value 1 /^*^ 
which the drop breaks is indicated by a star, corresponding to 
the turning point of the branch of solutions. Yellow line shows 
the interpolation of a discrete number of turning points. Inset: 
shape of the anchored drop at the critical value of 1/^*^. (b) 
Comparison between experimental values of 1/^*2 at which 
the drops break with the theoretical prediction from part (a). 
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We now turn to the breaking. The solutions to the 
elastica equation can be plotted as a family of curves in 
the (7^, 1/^^) plane, for different values of d, see Fig.[^. 
The curves all display a monotonic increase of Kq with 
1/^^, until a maximum value of where they reach a 
turning point, and after which no values of kq are found. 
The folding of the branch of solutions is associated to 
an exchange of stability at the marginally stable turning 
point. This implies that no stationary solutions to the 
elastica equation can be found beyond this point and 
the equilibrium can only be reached through a dynamic 
process which is not accounted for by the static model. 
We call the value of I at the turning point of the 
curve, and therefore expect that the droplet will break 
at Z = even though the calculated 2D droplet shapes 
the droplet neck still has a finite width, as shown in the 
inset of Fig. 

Experimentally, increasing quasi-statically the flow 
rate around a droplet pinned on an anchor amounts to 
walking along a curve with increasing This leads to 
an increase in kq until the value of ^ is reached, cor¬ 
responding to a critical velocity beyond which there 
are no stationary solutions. For each experiment we com¬ 
pare the fitted value of with the predicted value 
in Fig. [^. The agreement between the two is remark¬ 
able, confirming the interpretation of the breaking: above 
the velocity [/*, no equilibrium droplet shape can be 
found that satisfies the elastica equation with the im¬ 
posed boundary conditions. We find that the larger d/, 
the larger the value of 1/^*, i.e. in a given experimental 
setup with a fixed trap diameter d, the critical velocity 
t/* leading to break-up is smaller for large droplets than 
it is for small droplets. 

In many situations it is important to determine the 
volume of the trapped fluid in the anchor, for example if 
this is used to observe biological samples. This is equiv¬ 
alent to predicting, in our 2D view, the projected area 
Ad of the droplet left on the anchor after breakup. We 
expect the area to depend both on the diameter of the 
trap d and on the droplet radius R. The dependence of 
Ad on d is dictated by the the boundary condition that 
the trap imposes on the droplet shape. The dependence 
on R comes from the geometry of the break-up: we ex¬ 
perimentally and numerically observe that, for a given d, 
the pinch-off location increases with i?, occuring further 
downstream of the trap for larger drops. Rescaling all 
lengths by R, we show the evolution of the dimensionless 
area Adj^ as a function of d = d/R for both exper¬ 
imental results and theoretical simulations on Fig. 
Note that a perfect agreement cannot be expected be¬ 
tween the theoretical and experimental results since the 
model is 2D, whereas 3D effects are present close to the 
anchor. Still, the trends of the two curves are identical. 

We now examine the practical implications of the 
demonstrated method. A droplet breaks on an anchor 
as long as the outer flow velocity satisfies U > . In 


this sense, one does not need a precise flow control to 
produce a droplet: pushing the outer flow using a hand¬ 
held syringe is sufficient to break droplets. Also, the only 
relevant viscosity coming into play is the viscosity /i of 
the outer fluid: since does not depend on the inner 
fluid viscosity, the singularity occurs at the same time 
for droplets of a given volume, so that there is no dif¬ 
ference between breaking droplets of a viscous fluid such 
as FC-70 (/i = 24 cP) or breaking droplets of HFE-7500 
(/i = 0.77 cP). The volume left on the trap is dependent 
on the geometries of the trap and that of the drop: the 
sole important dimensionless parameter is d. Therefore, 
whatever the fluids used, their surface tensions, viscosi¬ 
ties, two drops of the same volume will break on an an¬ 
chor of diameter d into droplets of the same size, Ad^ that 
depends solely on d and R. 



FIG. 4. (a) Remaining area (dimensionless) as a function of 
d—djR. Dots: experiments, line: theoretical prediction. In¬ 
sets: two droplets of different sizes left on the same anchor 
(d = 150 /xm) by two different drops, (b) An array made up 
of 1568 individual drops produced through self-digitization. 
The aqueous solution contains fluorescein to aid in the visu¬ 
alisation. We have used square anchors with side d = 130 /xm 
in a channel of height /i = 35 /xm. Scale bar: 2 mm. 


The breakup of a single droplet can be generalized to 
an array of traps of any dimension, as shown in Fig. 
that shows an array of 112 x 14 drops. In this experi¬ 
ment an aqueous “puddle” initially fills the chamber and 
is then pushed by a flow of a wetting oil. The drops, 
which are produced at 10 Hz, detach as the front passes 
the anchors. Again, the interface connecting two anchors 
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is decribed by an elastica equation, albeit with different 
boundary conditions than the single anchor, and breaks 
when the enclosed volume is reduced below a critical 
value. 

This protocol can directly lead to applications to bio¬ 
assays, for instance by working with a suspension of cells 
in the aqueous phase. The droplet array recalls the mul- 
tiwell plate format, where each drop is analogous to a 
well encapsulating a cell population that is independent 
from its neighbors and can be continuously monitored. 
The oil phase can then serve to control the gas exchange 
with the drops, for example to oxygenate or de-oxygenate 
red blood cells, as shown previously [5]. 
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